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Abstract 

A new definition of the Wigner function for quantum fields coupled to curved space- 
time and an external Yang-Mills field is studied on the example of a scalar and a Dirac 
fields. The definition uses the formalism of the tangent bundles and is explicitly covariant 
and gauge invariant. Derivation of coUisionless quantum kinetic equations is carried out 
for both quantum fields by using the first order formalism of Duffin and Kemmer. The 
evolution of the Wigner function is governed by the quantum corrected Liouville-Vlasov 
equation supplemented by the generalized mass-shell constraint. The structure of the 
quantum corrections is perturbatively found in all adiabatic orders. The lowest order 
quantum-curvature corrections coincide with the ones found by Winter. 
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I. Introduction 



The Wigner-function method is weh-known as the most effective bridge between 
quantum and macroscopic regimes. It ahows one to derive from a quantum wave equation 
a Liouville equation for a quantum distribution function and explore the semiclassical Hmit, 
providing thus, a powerful ground for kinetic theory. Recently, the Wigner transformation in 
curved space-time was studied in different approaches0~@. In the present work, a covariant 
and gauge invariant Wigner function for quantum fields coupled to external gravitational 
and Yang-Mills fields is defined by using the horizontal lift of the derivative operator in the 
tangent bundle. The tangent and the cotangent bundles arise naturally in general relativistic 
kinetic theory and it was Vlasov, who first recognized thisi. This paper demonstrates, how 
the formalism of the bundles allows one to derive quantum Liouville equations explicitly in 
all orders of the perturbation theory. 



Sec. [I. recalls some features of general relativistic kinetic theory and gives the defi- 
nition of a covariant Wigner function, proceeding from considering expressions for dynamical 
observables (such as the current density vector, the stress-energy tensor and so on) of quan- 



tum fields. In Sec. |L| |AJ the first order formalism of Duffin and Kemmer is generalized to 
curved space-time and it is shown that a scalar field in the first order formalism is formally 
equivalent to a Dirac field. This allows one to derive quantum Liouville equations for different 



quantum fields in a universal fashion. Sec. III. derives quantum kinetic equations and gener- 



alized mass-shell constraints for the Wigner functions in terms of a semiclassical expansion. 
It is shown, by using the method of Ref. 0, how to eliminate the mass-shell constraint to get 



a transport equation for a reduced function on the classical mass-shell. In Sec. IV. some 
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unsolved problems are mentioned. The Appendices contain some technical details and useful 
formulae. 

Notations and conventions for a curvature tensor and a Yang-Mills field throughout 
are such that the commutator of two covariant and gauge invariant derivatives is 

' ~ ~ 1 Z6 

Vq,,V/3 Xfj_ = Raf3fiu X'^ + — Faf3 , (1.1) 

where e is a gauge charge and h is Planck's constant (the signature for a metric is (H )). 

II. Field equations and Wigner functions 

I will consider both a scalar quantum field and a Dirac field coupled to external 
gravitational and Yang-Mills fields. Let's consider first a scalar field. The model is described 
by the quantum field equation (see, for instance, Ref. 

(5°^(x)V,V^ + m^h^ - ^x)) ip{x) = VintM , (2.1) 

the Einstein equations: 

Ra(5{x) - ]^R{x)ga^{x) = k{T^^{x)) + kTl^\x) , (2.2) 

and the Yang-Mills equations: 

V<,F"^(x) = e{3^{x)) . (2.3) 

Here ^ is the nonminimal gravitational coupling constant, T^^\x) is the Yang-Mills field's 
stress-energy tensor and (J"(x)) and (Tq^(x)) are the ensemble averaged current density 
vector and stress-energy tensor of the scalar quantum field ip{x): 

(J«(^)) = (y)(v'^(^)^.^(^)), (2.4) 
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{T^pix)) = C-^f {^\x) v{x)) + 

+ (^^R^pix) + - 0(VaV^ - go.p{x)a)^ i^H^M^)) + (TL"i(x)) . (2.5) 

The last term in Eq. ( |2.5D depends on quantum interactions (the right hand side of Eq. ( |2.lD ), 
(p^{x) is a conjugate field and 

V„=V„ - , (2.6) 

the arrows indicate the directions of acting of the derivative operator. 

The ensemble averaging means that we have specified in-particles in an asymptotic 
past (with not necessary vanishing curvature, but admitting a reasonable definition of parti- 
cles) and a density matrix p based on in-states. Then the ensemble averaging of an operator 
A(x) is 

{A{x)) = Sp{pAix))/Spp. (2.7) 

Generally speaking, expressions like (|2.7D diverge and one has to renormalize it. It can 
be done in terms of a Wigner function by using, for example, a procedure like the BPHZ- 
procedure in quantum field theory0 or the dimensional regularizationil. We will return to 
this question elsewhere. In this paper all quantum operators throughout are considered to 
be somehow regularized. 

In the classical limit a system of scalar point particles is described by a distribution 
function F{x,p) which is a probability density in the phase space. The crucial feature of 
general relativistic kinetic theory is that in general relativity the structure of the phase space 
differs from the one in nongravity physics (see, for instance, Ref. ^). 

First, the second argument of a distribution function is a (co)tangent vector on a space-time 
manifold A4, which belongs to the (co)tangent space or the fibre T*{j\4) over point x (see, 
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for instance, Ref. 10). Coordinate transformations in the manifold, 

a;"'=x"'(x), (2. 



induce transformations in the fibre: 

Pa' = -g^Pa. (2.9) 

The set of cotangent spaces over ah points of the manifold forms the cotangent bundle over 
A4 with M. being a base space: 

T*{M) = U T:{M) , (2.10) 

and the distribution function F{x,p) is a scalar function on the base space and the cotangent 
bundle, that means that under the transformations ( |2.8D , (2.9) it behaves as the following: 



F{x,p) = F' [x'{x),^-p) . (2.11) 



An invariant volume element on the base space is known (see, for instance, Ref. 11), 

= (_^(x))l/2 ^3.0 ^ ^^1 ^ ^^2 ^ ^^3 _ ^2.12) 

Analogously, an invariant volume element on the cotangent bundle is 

dn(^) = (-c/(3;))-^/2 ^ ^ ^ ^p^ ^ (2.13) 

but the order of integration over the phase space in general relativity is crucial: to overcome 
an ambiguity, one should first integrate over the momentum space at fixed points x and 
then over the space-time manifold. In order to define procedures like coarse graining, which 
needs integration over finite space-time regions of a distribution function at fixed momentum 
variables p, one has to set a way of comparing the energy-momentum vectors of distant 
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particles. 

Secondly, the energy-momentum vector of a classical point particle, pa, is constrained by the 
mass-shell equation: 

5"^(x)p„p/j-m2 = , (2.14) 

where m is the particle's mass (we suppose that all particles are identical). More over, the 
energy po of a particle must be positive. So, the particle number in a phase space volume 
element is 

(IN = F{x,p)e{po)s{p'^ - m2)(in(^)dy(^) . (2.15) 

One can include 9- and (5-functions to the definition of the volume element on the physical 
sector of the phase space. Then after integration over po we will recover the 7-dimensional 
phase space. But as one will see below, for quantum particles the mass-shell constraint is 
more complicated and not only expressed by (^-function. It is thus more natural to define 
the classical distribution function as a distributed function, 

Mx,p) = F{x,p)e{po)5ip^ - m^) , (2.16) 

which has to be recovered from a quantum distribution function in the classical limit. If parti- 
cles couple to external (selfconsistent) gravitational and electromagnetic fields, the Liouville's 
equation which governs the dynamics of the distribution function is thenll 

(p"Z), - eF^p -^) Ux,p) = C[U] , (2.17) 

Opa 



D^ = Vc, + Tl^p^ — , (2.18) 



where 

d_ 

dpfs 

are the Christoffel symbols and C[f] is a collision integral depended on the mechanism of 
particles interactions. The operator Da is called the horizontal lift of the derivative operator 
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in the cotangent bundleS and it was introduced to kinetic theory by Vlasovi. It is invariant 
under the transformations 

The current density vector and the stress-energy tensor are expressed in terms of moments 
of the distribution function in the momentum space: 

-y=J=Pafcl{x,p) , (2.19) 

{T^p{x)) = I -^E= PaPMx,p) + {Ti$ix)) , (2.20) 



where the second term in Eq. ( 2.2C| ) depends on particles interactions. 



The goal of kinetic theory is to replace the field equation ( |2.lD by a kinetic equation 
for a quantum distribution function and express the expectation values ( |2.4D ,( |2.5D in terms of 
this function, so that in the classical limit (which has to be defined properly) we will restore all 
the set of Eqs. ( |2.16| )-( [2.2C| ). The most convenient method for exploring the classical limit is 
the Wigner-function one0. In the special relativistic quantum field theory the (oneparticle) 
Wigner function is the Fourier transform of a Green function with respect to the difference 
of its coordinatesEl: 

fix,p) = {TThy' J <fs e-2->'^/^(v^(xi)(^t(^2)) , (2.21) 

where 

^" = ^(^f-^f), (2.22) 
^" = ^(^? + ^?)- (2.23) 

In general relativity the situation becomes ambiguous because the definitions of the 
difference of coordinates ( p. 22 ) and a middle point ( 2.23| ) are not invariant under trans- 



formations of coordinates and there is no natural covariant generalization of the Fourier 
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transformation. The first attempt to generalize the definition ( |2.21| ) to general relativity 
and to evaluate the lowest order curvature corrections to the Vlasov equation was made by 
WinteJll. His definition is based on considering a geodesic connecting two points xi and 
X2- For evaluating the quantum corrections to the Vlasov equation in Winter's method one 
has to solve the geodesic and the geodesic deviation equations, which is a very complicated 
problem even in the lowest adiabatic order. 

A different approach! is based on introducing Riemann normal coordinates centered 
at any point in the vicinity of xi and X2, on defining the difference of coordinates and the 
middle point by Eqs. (2.22),( 2.23| ) in this coordinate system and on applying a certain dif- 



ferential operator on the Wigner function, in order to remove the dependence on the choice 
of a center of Riemann coordinates. This method is also covariant and is equivalent to the 
Winter's one, but it saves much work while deriving the quantum corrections to the Vlasov 
equation. 



An alternative approach was proposed in Ref. 14, which uses the formalism of the 
tangent bundles. This approach is explicitly covariant and also gauge invariant from the 
beginning and allows one to evaluate the quantum corrections to the kinetic equation up to 
any adiabatic order with few efforts. I will recall the main idea of Ref. 14. 

Consider first the expression for the current density vector ( [2.41 ). Let me rewrite this 
expression in the following way: 

(J„(x)) = ^ / A (^S\y)) {Tr^{x,-y)^^{x,y)) , (2.24) 

where is a tangent vector belonging to the tangent space 7^(7W) at a point x, 5^(y) is the 
4~dimensional (^-function, the trace is taken in the gauge group's representation space and 
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$(x, —y) and ^^{x,y) are the following series in the operators Vq, acting to the fields ip{x) 
and (f^x) at a point x: 

-y) = (^1 - y"V<, + VVaV;3 - Vv^V/jV^ + . . ip{x) (2.25) 

y)=(l + y'^Va + ^y°/v„V^ + Vv„V;3V^ + . • •) '^^(x) (2.26) 



The (5-function in Eq. (|]2|) can be represented by the standard way in a covariant manner, 
if one integrates over the cotangent space T*[M.) at a point x: 



Jt*{m) 



(2.27) 



Substituting the last expression into Eq. ( 2.24 ) gives the representation for the current density 
vector as an integral over the momentum space: 



(J„(x)) 



with the scalar function f{x,p) being 



V-9{x) 



Pa Trf{x,p) 



(2.28) 



f{x,p) = (TTh)-' J d'ye-^'y^-/" {^{x,-y)^^{x,y)) . (2.29) 



The expression ( 2.28 ) for the current density vector looks like the expression ( p. 19 ) in the 
kinetic theory, this indicates that we are on the right path. To be more sure, I will also 



represent the stress-energy tensor ( p.5|) in terms of the function (|2.29D . It can be done by 
the same way as for the current density vector by noting that an item 



(2.30) 



where the indices are symmetrized with — multiplier, can be represented by the same steps 
(|2l^)-(|2^) as the following: 

d'^ p 



V-9ix) 



PaiPa2 ■■■Pa^ Trf{x,p) , 



(2.31) 
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with the same function ( p. 29 ). After that, we can write the noninteracting part of the 



stress-energy tensor (^^) as an integral over the momentum space: 

(^afi{x)) = / , , . PaPp Trf{x,p) + 

+ + ( J - 0(V.V^ - 5«/3(x)a)) I -^E= Trf{x,p) + (T-j(x)) . (2.32) 

If one compares the expressions ( |2.28| ) and ( 2.32| ) with ( 2.19| ),( p.20| ), one can guess that in the 



classical limit the function Trf{x,p) will tend to the classical distribution function ( p. 16 ) 



The classical limit for the nonrelativistic Wigner function is usually explored from 



its transport equation (see, for instance, Ref. 15). So now, the next step is to evaluate the 



general relativistic quantum kinetic equation. But the definition ( |2.29| ) for a general rela- 
tivistic Wigner function is still useless because of the expansions ( p.25 ),( 2.26| ). The covariant 



derivative operator Vq doesn't annihilate the vector and one has to work a little harder 
to collect the series in Eqs. ( p.25| ),( p.26 ) into a convenient operator, as it was in nongrav- 



ity physics (compare with Ref. ^). To do that, let me introduce the horizontal lift of the 
covariant derivative on the tangent bundle0: 

V„ = V„ - r^^^ ^ (2.33) 



(compare with Eq. ( 2.18]) ). It can be easily shown that the operator does annihilate y 



./3 



and so, we can rewrite the n-th term of the series ( p.25D , (|2.26| ) in the following way: 



4 • • • • • • V„„ = ^ ( y"V, ) . (2.34) 



After that we get compact expressions: 



$(x, -y) = exp(-y°V«) (^(x) , (2.35) 
10 
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*t(a;,y) =exp(2/°V,)v9t(x) . (2.36) 
Eqs. (^),(|235D and ( pej) give us the definition of a covariant and gauge invariant Wigner 



function. As one will see below, this definition is very convenient for evaluating the quantum 
kinetic equation, but it looks rather formal if one wishes to compute the Wigner function 
explicitly for a concrete system. It is shown in [Appendix A that the Wigner function 



( 2.29| ) can be expressed in terms of base space quantities only, and conjunctions with other 



approaches are discussed. 

Let's now consider a Dirac field ip{x) coupled to external gravitational and Yang- 
Mills fields and satisfying the generalized Dirac equation (see, for instance, Ref. ^): 

(ih^'^ix) V„ - m) tP{x) = V.r^tm , (2.37) 

where 7"(x) are the Dirac matrices in curved space-time obeying the (anti)commutative rule: 

7"(x)7^(x) + 7^(x)7"(2;) = 25"^(x) , (2.38) 

and the acting of the derivative operator to a Dirac field is determined by the Fock-Ivanenko 
connection0r2cj(2;) : 

~ d ic 

V„ iIj{x) = — V'(x) - 0„(x)^(x) + -Aaix)^Pix) (2.39) 
ox°' n 

{Aa{x) is the Yang-Mills potential). 

A Dirac conjugate field ^^{x) = il^'^ {x)j{x), with 7(2;) commuting with all ^^"{x), satisfies the 
equation: 

i'ix) (^-in-f"{x) -mj =V~^], (2.40) 

with 

~ — d — — ic — 

V« i^{x) = — V(x) + V(x)rj„(x) - -V(x)A^(x) . (2.41) 

11 
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A covariant and gauge invariant Wigner function for a Dirac field can be defined by 
the same trick as for the scalar field. It is the following (matrix) function (the sign minus 
appears here because of the Fermi-statistics): 

N{x,p) = -{TTh)-' J d'ye-'^'y'^^-'^ (*(x,-y)*(x,y)) , (2.42) 

where the functions ^(x, —y) and ^(x, y) are related to the Dirac field tp{x) and the conjugate 
field ^/'(x) in the same way as for a scalar field (Eqs. ( |2.35| ),( p.36 ) and ( 2.33[ ) ): 



^{x, -y) = exp(-y°Va) il}{x) , 



"^{x^y) = exp(y"V„) V'(x) . 



(2.43) 



(2.44) 



In Ref. 14 the noninteracting part of the stress-energy tensor, the current density 
vector and the spin tensor of a Dirac field are expressed by using the generalized Gordon 



decomposition (see, for instance, Ref. 17), in terms of the independent components of the 



(2.45) 
(2.46) 
(2.47) 



Wigner function (2.42): 

A{x^p) = m^^ Tr N(x,p) , 

A''f^{x,p) = Um'^Tr {-i^'^{x)-t^\x)N{x,p)) , 

B{x,p) = —im~^ Tr N{x,p)^ , 

where 7^ = y/-g{x)eai3fi,.7°'ixh^{x)j^'{x)'y''{x). 

At the end of this section, I would like to note that a coupling to an external 
scalar field 0(x) can be described formally by changing the mass terms in Eqs. ( p.l[ ) and 



m I — > m + (j){x) . 



(2.48) 



12 
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Such coupling doesn't influence the definition of the Wigner functions but does infiuence the 
kinetic equation. 



A. The covariant DufRn— Kemmer formalism 

To derive the quantum kinetic equation it will be convenient for us to describe the 
scalar field within the first order formalism, in which the field equation looks formally like 
the Dirac equation. In special relativity this type of formalism was proposed by DuffinEl and 
developed by KemmeJli. I will extend it to general relativity. Let me introduce the new 
fields: 

Ua{x) = ihfi~^^'^Va^{x) ; U4(2;) = ^^/^(x) , (2.49) 



where fi is an arbitrary parameter of the mass dimension. Uniting the fields (2.49) into one 
five-component field: 



U4(x) 



and defining the 5x5 -matrices: 



(2.50) 



r"(a;) 



O 



[g^f^ix] I J 



(2.51) 



with 



M{x) 



o 



o 



m{x) J 



m{x) = [m — $^fi R{x)) , 



(2.52) 



(2.53) 
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one can rewrite the equation (^]^) in a form which is formany the same as for a Dirac field 
( 2.37] ) (with a matrix variable mass): 



(inr^ix) v„ - M{x)) u(x) = y/„ju] , 



(2.54) 



the only difference is that the matrices r"(x) obey the Duffin-Kemmer commutative rule 
rather than the Dirac one ( |2.38| ): 



r"(x)r (x)r'5(x) + r'^(x)r'^(x)r"(x) = g''''{x)T''{x) + g'"'{x)T''{x) . (2.55) 



The covariant differentiation in Eq. ( 2.54 ) is defined by the connection Aq,(x 



K{x) 



SO that (compare with Eq. ( p. 391 )) 



Vq, u(x) = u(x) - Aq,(x)u(x) + —Aa{x)\y{x) 



h 



If we define the matrix r(x), which commutes with all r^°(x): 



r(x) 



^ g^f^{x) I O ^ 



O 



then the conjugate field 



u(x) = u^(x)r(x) = (u^"(x),ui(a 



satisfies the following equation: 



(2.56) 



(2.57) 



(2.58) 



(2.59) 



u(x) -ihT^{x) V„ -Mix) = Vl^M 



(2.60) 



14 
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It can be checked that the covariant derivative operator annihilates the F-matrices (p.51| 



and ( 2.5^ ) as well as the Dirac 7-matrices: 



d 



v^Tf'{x) = ^T^{x)- A„(x),r^(x) +ri(x)r(x) = o, 



d 



V, T{x) = — T{x) - Ki{x) T{x) - r(x) A°(x) = . 



(2.61) 



(2.62) 



The analogy between the scalar field in the first order formalism and the Dirac field will 
show itself more explicitly if one looks at the expressions for the commutator of two covariant 
derivatives acting to each field. If one denotes: 



Us(x) 



u(x) for s = 1 (the scalar field) 
^{x) for s = 1/2 (the Dirac field) 



and 



r?(x) 



r"(x) for s = 1 

7"(x) for s = 1/2 
then it can be obtained easily from Eqs. ( p.39D , (|2.57| ) that 



(2.63) 



(2.64) 



Us(x) = Asapix) \is{x) , 



(2.65) 



with 



le 



For the conjugate fields one gets: 



(2.66) 



Va,V/3 \ls{x) = -\ls{x) Asaisix) 



(2.67) 



Let me now define the unified Wigner function: 



Ns{x,v) = es{Trnr^4-9^) / d'ye~2*^"P"/"(U,(x,-y)U,(x,y)) 



(2.68) 
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with \Js{x,—y) and \Js{x,y) being related to Us(x) and Us(x) respectively by means of 
formulae analogous to (^351), (|23^) and 



(2.69) 



1 for s = 1 
-1 for s = 1/2 

Then for s = 1/2 one will get the Wigner function of the Dirac field ( 2.42| ) immediately. 
In order to extract the Wigner function of the scalar field ( |2.29| ) from ( 2.68 ), one should 
introduce the projection operators: 



Pi = ^ (r-(x)r„(x)-i) 



7^2 = 3 (4-r-(x)r,(x)) 



o 



o 



5P 



o 



\ 



(2.70) 



1 ; 



(2.71) 



with the following properties: 



Then 



Vi+V2 = l ; Pi 7^2 = ; Vt=V^ 



Vi r"(x) = r°(x) V2 ; V2 r"(x) = r"(x) Vi , 



r"(x) r^(x) Vi = c/"^(x) Vi . 



fix,p)=^i-Hr {ViNiix,p)) 



(2.72) 

(2.73) 
(2.74) 

(2.75) 



In conclusion, it is worth noting that the Duffin-Kemmer formalism can be intro- 
duced for a vector field also. The evaluation of a quantum kinetic equation in this case was 
carried out in Ref. |^. 
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III. Quantum kinetic equations and 



mass— shell constraints 

A. The first order equation for the unified Wigner function 

If one takes into account the formal analogy between the Dirac field and the scalar 
field in the first order formalism, one can guess that the equations which govern the dynam- 
ics of the Wigner functions, look formally similar for both fields. These equations can be 



extracted from the equation for the unified Wigner function (2.68). Before proceeding to the 



derivation of this equation, I will prove the following identity (see also Ref. pO|): 



De,Ns{x,p)=es (7rn)-y-g(x) j e-"'^^^"^/^ V«(U,(x, -y)U,(rE, y)) , (3.1) 
where Da is the horizontal lift of the derivative operator in the cotangent bundle ( Eq. ( p. 18 ) 



with the minimal gauge invariant extension ) and Vq, is the horizontal lift of the derivative 
operator in the tangent bundle ( 2.33| ). 



The identity can be easily proved if one works in the Riemann normal coordinates 

centered at point x0. In this coordinate system both the operator Da and the operator 
Vq coincide with the gauge invariant extension of the partial derivative operator da and the 
identity ( |3.1| ) is trivial. The proof will be completed if one takes into account the invariance 



of the derivative operators D^ and and of the volume element in Eq. (3.1). 

Utilizing now the identities (|B2[) , (p^) from Appendix B , we will get, after integrating once 

by parts: 

2i 

DaNs{x,p) = —paNs{x,p) + 

+ 6, {nhr^ ^-g{x) Jd'^y e-^'^^^/^ (2 (^e^^"^- V„ u,(x)) U,(x, y) - Gsa{x, y)) , (3.2) 

17 
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where 



Gsa{x, y) = Us(x, -y) (Caix, y)Usix, y)^ 
(^^(x, -y)\Js{x, -y)) U.(x, y) + 2 (h^{x, -y)U,(x, -y)) tj,(x, y) . (3.3) 



So far, the equahty (|3.2| ) is an identity. In order to obtain an equation, let me 
multiply it by iTi/2T^{x) ( with a summation over a ). Taking into account the property 
( 2.61| ) and using the field equations ( p. 37 ) or ( |2.54| ) together with the identity (B5), one will 



get: 



M(x, -y ^) Ns{x,p) = r-(x) [p^ + '^b^ ) Ns{x,p) + 



+ y C(x) e, {TTfir^ ^-g{x) jd^y e-'^^'^v^l^' (G,,(x, y)) + Cf , (3.4) 
where M{x, — given by Eq. (|B^ ) with the substitution: 



and Cf°' is a "collision" term which depends on the interaction terms in Eqs. ( 2.37] ) and 



( 2.54| ). I will not discuss this term in this paper but return to this question elsewhere. In 



the rest of the paper I will omit the "collision" term and deal only with collisionless kinetic 
equations. 

The equation ( |3.4| ) will be a kinetic equation if one also expresses the second term on the 
right-hand side in terms of the function Ns{x,p). It can be done by using the results of 
Appendix C . Substituting Eq. ( |C12| ) into the integrand in Eq. ( |3.4D , using the identity (3.1) 



once again and integrating by parts, one will get the following (collisionless) equation: 

M{x,-'-^-^)Ns{x,p)=r^{x) (^^LD. + 4^„p.)iV,(x,p) + 

+ yr^(x)i,i, N,{x,p) + yr^(x) Ns{x,p) As2a ■ (3.6) 
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Here -Ri 2« ^si,2o are represented by the infinite series in the momentum derivatives: 

oo 

Rho. = C + Y: KlfRf^' , (3.7) 

k=2 

oo 

M2. = T.KlfA%-\ (3.8) 

k=2 

with 



oo oo 



kZ = 5f CI + E E • • • E Ch,...,Ji'"^^" . . . R^^f (3.9) 

n=l fci=2 k„=2 

and C^^ being numerical coefficients (their exphcit expressions are presented in [AppendEi 
C). The operators Ra ''^ and Aia in Eqs. (|3.7|)-(|3^ are defined as follows: 



^ (^)fc <,,,^.3...,, aV9p.. ...Pu,, (3.10) 
^i'^ = (y)'^a.i;....... d'^/dp., ...Pu,, (3.11) 



with being defined by ( 2.66 ) 



The equation (|3.6|) together with the definitions (|3.7|) - (|3.11|) give us the quantum 



kinetic equation for the unified Wigner function (|2.68| ) in the collisionless approximation. But 
this equation is not yet in a useful form for exploring quantum kinetic properties. Firstly, 



the infinite series in Eqs. (3.7)- (|3.9| ) are formal as long as nothing is known about their 
convergence. The situation here is quite similar to the one arisen in the nonrelativistic case 
(see, for instance, Ref. |15|). One has to assume that quantum interference is suppressed 
somehow and that the external fields are "slowly varying" at the effective Compton scale 
(compare also with the discussion in Ref. |2|). Then Eqs. (|3.7])-( |3^ can be approximated by 
the first few terms and Eq. ( |3.6[ ) will be a differential equation of a finite order. On the other 
hand, in some cases the infinite series in the momentum derivatives can be collected into 
integrals and the kinetic equation will be an integro-differential equationEl. If the quantum 
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system is in a highly coherent state then some selective summation of infinite subseries could 
provide a reasonable approximation^. The problem is a serious one but its discussion will 
be reserved for further study (for the recent discussion see Ref. |2^). It is worth noting that 
similar expansions arise in the description of the extended bodies in general relativityil. 
Secondly, the equation (^) doesn't look like a transport equation. It will do so if we " square" 



it (compare with the evaluation of the transport equation for 1/2-spin particles in Ref. 13). 



For this purpose let me write Eq. (3^) in compact notations: 



MN,{x,p) = K,N,{x,p) , (3.12) 

where by Kg I have denoted the operator acting on the function Ns{x,p) on the right-hand 
side of Eq. 

Later on, the scalar field and the Dirac field will be treated separately. 
B. The scalar field case 



An equation for the Wigner function of the scalar field ( |2.29 ) can be obtained by 



the means of the prescription ( 2.75| ). Multiplying the equation ( p. 12 ) from the left by the 



projection operators ( 2.70| ) and ( ^.71 ), one gets, by using the properties ( 2.73 ) and the explicit 



matrix form of Eq. (p.q), the following two equations: 

ViMNi{x,p) = KiV2Ni{x,p) , (3.13) 

V2MNi{x,p) = KiViNi{x,p) . (3.14) 

Taking the structures of the matrices ( |2.52| ) and ( 2.70| ),( p.71| ) into account, one can easily 
obtain 

Vi M = 11-^ mVi , (3.15) 
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V2M = ijV2, (3.16) 

where the operator m is given by 



Using ( p.l5| ) and ( p. 16 ) in ( p.l3| ) and ( p. 14 ) leads to the following equation: 



- ^h''R{x, - ri Ni{x,p) = Ki Ki Vi Ni{x,p) . (3.18) 

Substituting the explicit form of the operator Ki from Eq. (^|^) into Eq. (|3.18|) and taking of 
the property ( |2.74| ) into account, one can obtain the equation for the Wigner function ( p. 291) 



of the scalar field by making use of the prescription ( p. 751) . The result is: 



m' - ih'R{x, - j f{x,p) = [K'' + -Ri^ j f{x,p) , (3.19) 

where is the following operator: 

Ka f{x,p) = (^R-.^D, + R!^^p, - ^Fi^) f{x,p)-^f{x,p) F2a, (3.20) 



and the operators Ria and Fi_2a are defined by Eqs. ( |3. 81 )^( ^1X1 ), replacing: 



Rau for Rla 

Asau ^ { ■ (3.21) 

Fau for Fi,2a 

Eq. ( p.l9| ) contains, in fact, two equations. To see that, one should remember that the Wigner 
function is an hermitian (matrix) function in view of the definition ( 2.2S| ). Separating the 



hermitian and the antihermitian parts of Eq. ( |3.19| ) , one will get two independent equations 
for the only function f{x,p), which have to be satisfied simultaneously. Here I will only 
write out the explicit form of the classical terms and the first quantum corrections of these 
equations. 
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The first quantum corrections due to the Yang-Mills field are multiplied by h while the ones 
due to the gravitational field are multiplied by h"^. It is worth noting that to get the terms of 
the order of in both equations, one has to expand Eq. ( |3.19D up to the terms of the order of 
^n+i^ because the antihermitian part of it includes the overall multiplyer ih (we suppose that 
quantum interference effects have been quenched by averaging, so that the Wigner function 
itself is independent on h, see also the discussion after Eq. (|3.1l| )^). Then two equations to 
the single function will be consistent up to the terms of the order of h"' (see the discussion 
below) . 

These two equations are 

~ P f) 1 Pi) (J ' ~ 

p-D^f{x,p) + -p''—{F^,,f{x,p)} = — — \F,^,D''f{x,p) 

2 OPa 8 Opa L 

The [, ]-bracket stand here for the commutator and the {, }-bracket - for the anticommutator 
of two matrices in the gauge group's representation space. Remember also that the operator 
Da is defined as following: 

D^f{x,p)= [da + Tl^p^^] f{x,p) + -[AaJ{x,p)] . (3.24) 

Eqs. (3.22) and ( |3.23D were first obtained by WinterHfl through a very different method. The 



first equation gives the quantum corrections to the mass-shell constraint of classical parti- 
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cles ( p. 141) and the second equation is the quantum Liouville-Vlasov equation. These two 
equations must be satisfied simultaneously. It could lead to inconsistency with the classical 
limit if we were not able to reduce them to a single equation. 

Indeed, in the naive classical limit (h — > 0) the Wigner function satisfies two classical equa- 
tions: 

Qf{x,p)=0, (3.25) 
Cfix,p) = 0, (3.26) 

where 

n = p^p^ - (3.27) 

and L is the classical Liouville operator, which acts on the left hand side of Eq. ( |3.23| ) ( it 
coincides with the left hand side of Eq. ( 2.17| ) if the gauge group is U{\) )• One can easily 



check that the operator L annihilates $7 and therefore the function 

fd{x,p) = F,i{x,p)b{^) (3.28) 



will satisfy Eq. ( p. 26 ) if Fci{x,p) satisfies the same equation but on the mass-shell only. It 



means that the following equation holds: 

£.Fa{x,p) = nA{x,p) , (3.29) 

where A{x,p) is an arbitrary function, which is nonsingular on the mass-shell. 

We see that in the classical limit one is able to reduce two equations for the Wigner 
function (which is a distributed function in the relativistic case) to a single equation which 
coincides with the Vlasov equation for the 7-dimensional distribution function on the mass- 
shell ( more exactly, the sum of the distribution functions for particles and antiparticles, 
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compare with Ref. 13 ). The same situation can be expected in each order in h, because 
otherwise one wouldn't get the right semiclassical hmit. 

To show that Eqs. ([3 .221) and (|3.23|) lead to a single equation up to the terms of the next 
adiabatic order, I will use the method of Ref. |2|. Let me look for a solution to Eg. ( 3.22| ) of 
the form: 

f{x,v) = Fo{x,p) 6{Q) + 6'{Q) + F2{x,p) 6"{Q) , (3.30) 

where Fq{x,p) is the quantum corrected distribution function on the mass-shell: 

Fo{x,p) = Felix, p) + Fqu{x,p) , (3.31) 

and Fi{x,p) and F2{x,p) have purely quantum origin. It is worth noting that if one has been 
interested in all terms in the adiabatic expansion then, as it is mentioned in Ref. ^, one should 
include the term ?i(^ — l/6)R{x) in to get an improved asymptotic approximation for the 
Wigner function. But as, for the moment, we are only interested in the lowest adiabatic 
order, it is more convenient to consider that this term also contributes to the ojf mass-shell 
part of the Wigner function. 
Using now the property: 

fi5(")(0) = -n5("-^)(17) , (3.32) 
and treating all derivatives of the (5-function as linearly independent, one can get easily 



from Eq. (3.22) the first quantum corrections to the classical distribution function ojf the 
mass-shell: 

Fi{x,p) = UF,i{x,p) , (3.33) 



1 r - 



h 



2 



F2{x,p) = -U,n F,i{x,p) = -—Rf,,p^'p''F,i{x,p), (3.34) 



3 
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where 11 is the operator acting on the right hand side of Eq. ( p.22| ) . 

The substitution of the trial solution (|33c| ) with (|33l|) , ( |333D and (|]3|) into Eq. (|]23|) gives, 
after omitting the higher adiabatic order terms, 



6{n) C (Felix, p) + Fqu{x,p)) + 6'{n) Cfl Felix, p) + -6"in) t 



6in)AF,iix,p)+6'iQ) 



Fciix,p) + -5"in) 



A,0 



Felix, p) = 

Felix, p) , (3.35) 



where A is the operator acting on the right hand side of Eq. ( 3.23| ) (the operators ft and A 



are written out explicitly in Appendix D ). 

Eq. ( p.35| ) seems to give rise to inconsistency of the method, but hopefully, the structures of 
the operators in Eqs. ( |3.22| ) and ( 3.23| ) are so fine tuned that it leads to the cancelling of the 



terms which multiply the derivatives of 5-function. Indeed, as it is shown in Appendix D, 
the operators ft and A satisfy the following identity: 



Using also the fact that the operator C annihilates 17 and, therefore. 



(3.36) 



(3.37) 



one can transform Eq. ( 3.35| ) into: 



5in) t iFelix,p) + Fe^uix,p)) + b'i^) ft£ Felix, p) + 



t Felix, p)=5ia)k Felix, p) 



(3.38) 



Remember now that the classical distribution function Feiix,p) satisfies Eq. ( p.29| ). This 
leads, after using the property ( p.32| ), to the equation for evaluating the first quantum cor- 
rections to the classical distribution function on the mass-shell: 



5in)tFguix,p) = 5in) (kFeiix,p) + f[^ix,p) 



(3.39) 
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Eqs. ( |3.3C1| ),( 3.31 ),( 3.3^ ),( |3.34| ) and ( p.39|) solve, in principle, the task of finding the lowest 



order quantum corrections to the classical distribution function ( 3.28| ) of scalar particles due 
to the coupling to external gravitational and Yang-Mills fields. As it was emphasized in Ref. 
^, some terms of the higher order must be added, they give the right trace anomaly of the 
stress-energy tensor. 



C. The Dirac field case 



An equation for the Wigner function of the Dirac field ( 2.42 ) can be obtained by 
squaring Eq. ( 3.12| ) for s = 1/2. Remembering that in this case the operator M coincides 
with the field mass m, one easily obtains: 



m^N{x,p)=ki/2Ki/2N{x,p) . 



(3.40) 



Utilizing the algebra of the Dirac matrices ( 2.3§| ), one can deduce from Eq. (3.40) the equa- 
tions for the independent components of the matrix Wigner function It happens 
that the equations form a coupled system of equations and other components of the Wigner 
function ( |2.42| ) ( a vector and a pseudovector functions ) defined by: 



A'^{x,p) =Tr{-f''N{x,p)) , 



(3.41) 



5"(x,p) =rr (7V^(a;,p)) , (3.42) 

are expressed in terms of the functions ( |2. 451 )^( ^.47 ) and their derivatives only. 
The explicit form of the equations is very unwieldy in general and here I will write out the 
equations with the first quantum corrections only ( of the order of /i ). The general form of 
the higher terms in expansions like ( p. 19 ) for the Dirac field can be found in Ref. |^. It is 
only worth noting that quantum corrections of two types arise. The corrections of the first 
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type are similar to and coincide, up to the terms of the order of with the ones in Eq. 
( 3.22| ),( ^l23| ) when ^ = 1/4. The corrections of the second type reflect the coupling of the 
spin and/or the isospin of quantum particles to external gravitational and Yang-Mills fields 
and they mix different components of the Wigner function (^)~ (pl7|) . 



As in the scalar field case, Eq. (|3.40D leads to two systems of equations: the constraint 
equations, which describe the quantum corrections to the classical mass-shell constraint 



( 2.14 ), and the quantum transport equations governing the evolution of the Wigner function. 
These equations with the lowest order quantum corrections are 



{n + U)A{x,p) = f{F^,,A''''{x,p)} 

{n + U)A^,{x,p) = f {F^,,^(x,p)} ---^{f^,,B{x,p)} , 



(3.43) 



(3.44) 



{C - A) A{x,p) = -S^, A^^-{x,p) 
{C-A)B{x,p) = -§^,A^''^{x,p) 
{C - A) A^„{x,p) = -\S^y A{x,p) + \S^y B{x,p) , 
where and the operators 11, C and A are defined as in the previous section (they are written 
out in Appendix D| , only the terms of the order of 7i of 11 and C are kept in Eqs. ( |3.43| ), (|3.44 ) 
) and Sfj_i, is the following operator: 

S^^A{x,p) = - [Ff,^,A{x,p)] + —Rp^^y—A{x,p) + — — ^(x,p)} . (3.45) 

The tilde over tensors indicates dual tensors, for example, F^^ = ^y^—g{x) e^uap F"^ . 

Eqs. ( p.43D ~( ^^ ) were obtained for the case of the C/(l)-group ( an external elec- 
tromagnetic field ) in Ref. |l^. It is interesting to compare these equations with the heuristic 
equation proposed by Israelii for describing particles with internal spin in external gravi- 
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tational and electromagnetic fields. He introduces the extended phase space by adding to 
it the spin angular momentum of particles, and the distribution function N{x,p,s) is a 
function on this space. One can show that the transport equation to the function N{x,p, s) 



leads, after integrations over s^^'^ , to the equations ( 3.43 )~( 3.45 ), if one identifies the functions 
A{x,p), B{x,p) and A^'^{x,p) with the corresponding moments of the function N{x,p,s) in 
the s^'^-space. The correspondence here is quite similar to the one which arose in the case 
of colored particles in a Yang-Mills field ( see, for instance, Ref. ). 

IV. Final remarks 

As a conclusion, I would like to mention a few further developments related to 
the present work, which have not been considered in the paper. (1) This paper didn't 
specify any regularization procedure as well as it didn't discuss in detail the question of the 
renormalization of the statistical averaged operators. This question is especially important 
when one wishes to describe the particles creation due to quantum interactions or external 
fields. This process could be described in kinetic theory, possibly by a source-term in the 
Liouville equation and could lead to Boltzmann entropy generation ( different approaches 
to defining a nonequilibrium entropy for quantum fields in a cosmological setting have been 



proposed in Refs. 28 -3C, see also Ref. E). (2) Specification of particles interactions will 



give an explicit expression for the collision term in Eq. ( |3.4| ). As it is shown in Ref. 31, a 
satisfactory approximation is provided only by the two-loop order of the perturbation theory 
and higher. The theory becomes very involved in this approximation and one needs some 



physical hypothesis to make further progress. (3) In Ref. 32 the quantum kinetic equation 
for the Wigner function of a scalar field has been solved explicitly in the lowest adiabatic 
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order in a Robertson- Walker space-time, in Ref. ^ - in a class of manifolds admitting a 



Killing vector ( for instance, static space-times ) and in Ref. 34 - a conformal Killing vector 
( the later class includes, in particular, a Robertson- Walker space-times ). Ref. ^ considers 
also a few other physical systems permitting an explicit analysis. In each case the quantum- 
curvature corrections to the local-equilibrium distribution are expressed by local geometrical 
quantities and a few momentum derivatives of the classical distribution function. This fact 
extremely simplifies the Einstein- Vlasov problem, reducing it to an analysis of differential 
equations of a finite order ( such analysis by using the moments method was carried out for 



the Friedmann cosmology in Ref. 35). A natural question arises: is this property general or 
is it only a lucky chance and how do space-time symmetries show themselves in the Wigner 
function? The better understanding of the structure of the relativistic phase space could 
possibly clarify this tie. 
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Appendix A The Wigner function as expressed on the base 



space 



The definition of a covariant and gauge invariant Wigner function ( p. 291) is not very 
useful if one wishes to compute it exphcitly or explore its ultraviolet behavior, because the 
fields y) and ^^{x,y) have been defined on the tangent space rather than on the 

physical space-time. Fortunately, one can express these fields in terms of fields on the base 
space, though in general the expressions also include infinite series in the derivatives. 
Consider first the field —y)- Let me rewrite expression ( |2.35| ) in the following form: 



= exp(-y°Vc,) exp{y°'da)) iexp{-y°'da) (p{x)) 



(Al) 



The second cofactor in ( [Al| ) is obviously ip{x — y) for any numerous value of in each fixed 
coordinate system. The first cofactor is in general a differential operator acting to ip{x — y). 
It can be written more explicitly by using the Campbell-Baker-Hausdorff formula ( see, for 
instance, Ref. |3^ ): 

exp(i) exp(B) = exp , (A2) 



\n=l 



with Cn being found by the sequence: 



Ci = A + B 



C2 = - 

2 



A,B 

n-l 



Cn — — :~ \ -A, 



1 1 



n ' — '„ ml 

m=2 



A, B 



fci 



Cfc2 ) C'fcl 



(A3) 
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If one puts first, B = y'^Va and A = —y'^Va = —B — ^y°'Aa{x), then it follows from Eqs. 



that 

Z{x, -y) = exp(-y"Va) exp(y°Va) = 

(OO / 2 \ 

This expression can be written in a compact form by using Feynman's path ordering operator 
P: 

-v) = P exp (^-| f; ^y-" • • • y-"^„o..,...„„ j . (A5) 

In the same way, with B = y°(9„ and A = -y"V„ = -B + Tl^„y°'y'' ^ Eqs. (^,(^ give: 



'Pix, -y) = exp(-2/°V„) exp(y"aa 



= exp [y-y^K, [^-d. + — j - y^y^^K,^, + 2 + " ' 'J ' ^^'^ 

If one uses the memann normal coordinates^ centered at point (RNC), the operator 
V{x, —y) vanishes. Then the field —y) differs from (p{x — y) only by a phase, which is 



In the same way, the field in RNC equals to {x,y)(l)^ {x + y) and the Wigner 



function (^.29) becomes 



f{x,p)\RNC = (Trh)-' / d'ye-''y''P-/^Zix,-ymx-y)cl,\x + y))zHx,y) (A7) 



In absence of an external Yang-Mills field, Z{x, —y) = 1 and the function (A7) coincides 
with that defined in Ref. § ( in RNC ) within a factor equal to the square root of the Van 
Vleck-Morett determinant, or the one defined in Ref. |l|. In other coordinate systems the 



Wigner function can be calculated by using Eq. (A6) ( see also the discussion in Ref. ^ ). 
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In another limiting case, when the space-time is flat, the tangent space coincides with 



Minkovski space. In this case the function (A5) takes the familiar form: 



Z{x, -y) = P exp 



dty'^Ao,{x-ty) 



(A8) 



and the gauge-invariant Wigner function (^3) coincides with that defined in Ref. 25 



or m 



Ref. 27, in the Dirac field case (for the case of the gauge group being ^7(1), see for instance 



Ref. 

Thus, our definition of the Wigner function in curved space-time is locally equivalent 
to those of Refs. But for nontrivial topologies, the different definitions of the Wigner 
function for a system in a highly coherent state could lead to globally different results. 
The same analysis can be carried out for the Wigner function of a Dirac field ( |2.42| ). 



Appendix B Useful identities 



In this Appendix I will prove a few useful identities, which have been used in deriving 



the equation for the unified Wigner function (3.6). 
Let A and B be two arbitrary operators. Then the following well known identity holds: 



B 



E 

n=l 



1 



B , A 



Putting A = Vq and B = —y^Vy leads to 



(Bl) 



V, U(:e, -y) = e-^"^" u(x) - i?„(x, -y)U(x, -y) , 



(B2) 



with 



H,{x,-y) = J2^—j^y''K..y''- 



n=l 



nl 



V 



1^1 ; 



(B3) 
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The field U(x, —y) is defined by the same formula as in (|2.3 



U(2;, -y) = exp(-y"V„) u{x) . 



(B4) 



If one puts A = M{x), then (Bl) with the same 13 gives the following identity: 



g-jz-v. M{x)u{x) = M{x, -y)lJ{x, -y) , 



(B5) 



with 



CO f_-i\n 

M{x, -y) = Mix) + J2 -^y"' • • • y"" V,i • • • V,„M( 



X) 



n=l 



(B6) 



The next identity follows from the formula for the left derivative of the exponent of an 
operator function B{t): 
/ 



d 
— f 

dt 



V 



dB{t) 
dt 



1 



■ U(x, -y) = -Va 'U{x, -y) - G„(x, -y)U(x, -y) , 



where 



(-ir 



- (n + 1)! 

The same identity for the conjugate field is obtained by reversing the sign of y° 



(B7) 



Substituting B{y°') = — i/'^Vq, — Ylu^taV'^'^ f (with no summation over a) into ( pTI) gives the 
identity: 



(B8) 



(B9) 



Appendix C Evaluation of Gsa{x,y) 



In this Appendix I evaluate expression with the operators H and G being 

defined in (|B3D , ( p9|) . I will omit here the subscript s for convenience. 
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Let me first prove the following recursive formula (n > 2): 



d 



n-2 



4") = -4-1) U(x, -y) + R^:> — U(x, -y) - ^ C,"" ^ R{n~k)u ^(fc) (^1) 



where = n\ /(n — h)\ k\ are the binomial coefficients and 



=y--... y^- 



U(x, — y) if n > 1 



(C2) 



Gi') = V«U(x,-y) , 



R^S^^ = y''^ ...y^-Rlr.u.-u. if n > 3 



A. 



(n) 



. . . y^" A 



if n > 2 



-^a — y -^o 



(C3) 



(C4) 



with Aau being defined by Eq. (2.66). 



For n = 2, the last term in Eq. (CI) disappears and the validity of this formula follows from 
the definition ( 2.33| ) and Eq. ( 2.65| ). Suppose that Eq. (CI) holds up to some fixed n. Then, 
for n + 1, one can get the equality: 

= -4") U(x, -y) + 4"+!)^ A u(x, -y) - 

n-2 

E' 

k=l 



- R^^> U(x, -y) - J2 Cfc"^ (4""''^^^'' Gi'^^ + 4"-'=)'^ Gi^'^+i)) . (C5) 

A:=l 

The proof will be completed if one takes into account the properties of the binomial coeffi- 



cients: 



C" = l. 



(C6) 



Let me now define the series: 



oo 



n=2 



(C7) 
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with arbitrary ^f^'s depending on n. A repeated use of the formula ( |Cl| ) gives: 



+ E E • • • E (-1)" c!r' ■ ■ ■ cz:^ Rt^^- ■■■Hp, 

n=2fei=2 fc„=2 



X - + (^1 - mtr^'^,) U(x, -y) , (C8) 



where 



Sm = ki + . . . + km ■ (C9) 

By putting gn = [—l)'^/n\ and gn = (— l)"/(n + 1)! , one wiU get expressions for the acting 
of the operators ( p3|) and (^9| ) to U(x, — y) and then evaluate the function Ga{x,y), (^). 
Unfortunately, the expression for Ga{x,y) includes terms with 

(af° U(2;,y) or ^Vq,U(x, — y)^ tj(x,y). In order to obtain an expression includ- 

ing only operators acting to U(x, —y)Xj{x,y) , one should use the identity ( |B8| ) for field U 
as well as for the conjugate field. It is easy to get: 

^ ^ ■ U(x, -y)) U(x, y) = l^ (U(x, -y)U(x, y)) - (U(x, -y)U(x, y)) - 



- iu(x, -y) y)U(x, y)) - ^ -y)U(x, -y)U(x, y)) , (CIO) 

(^V„ U(x, -y)) U(x, y) = (U(x, -y)U(x, y)) + (U(x, -y)U(x, y)) + 

+ ^U(x, -y) y)U(x, y)) - ^ -y)U(x, -y)U(x, y)) . (Cll) 

A repeated use of the identity ( |C8D gives then, after elaborate calculations, the formal rep- 
resentation for the function Ga{x,y), (|3.3| ) as an infinite series in operators acting to the 
product UtJ = U(x, — y)U(x,y): 

Ga{x,y) = (^i^r,^. - K^a-^ + UU + UU^2a , (C12) 
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where 

-1,2/3 Ty{k)u 



k=2 



fc=2 



with 

oo c» oo 



= '^f Q + E E • • • E ^^L....^^'^"^^" • • • 4:^' , (CIS) 



n=l fci=2 A;„=2 

and da^'^ and ^i*'^ being defined by (|C|),Q. 

The numerical coefficients C|.^ can exphcitly be found in the fohowing way. Let A{t) be 
an operator function. Let's introduce the fohowing generating functionals: 

oo oo oo 

^i^] = E E • • • E cL..k.A^''^ ■ ■ ■ ^^'"^ , (C16) 

n=lfci=2 kn=2 

where A'^^'> = d}'A{<d)/dt^ are the derivatives of A{t) at t = 0. 

Then, the coefficients C\_^ are found if we define these functionals as follows: 

K^[A] = -l + KB, (C17) 

K\A] = -l + KC, (C18) 

K\A] = -2C' C-^+KC C-^ (C-l) , (C19) 

K^[A] = -K (C-l) , (C20) 

where 

K = (l- B' B-^ + C'C-^) (^^BB-^ + ^CC-^^ . (C21) 

Here B = 5(1), C = C(l), B = B{-1) and C = C(-l), where B{t) and C{t) are the 
following functionals: 

oo oo oo n -, 

m = ' + E E ■ . . E (-1)"''" n . («^) 

n,=lfci=2 fc„=2 1=1 ^ * / IV J I / 
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OO OO OO IL -I 

B' and C" denote the derivatives at t = — 1. 

The functionals B{t) and C{t) can be written in a compact form by using Feynman's path 
ordering operator P: 

B{t) = P exp (^J^dxb{x)^ , (C24) 
C{t) = P exp (^J^dxcix)^ , (C25) 

with 

b{x) =x-^ r dzz^A^^^\z) , (C26) 







JO 

In particular, the coefficients in Eq. ( C15| ) needed only for evaluating the linear terms with 



respect to the curvature tensor in Eq. ( |3.6| ) are 

d = (-l)Vfc! + (l - (-1)^) /{k + 1)! , (C28) 

C| = (-l)Vfc! - (l + (-1)'=) /{k + 1)! , (C29) 

C| = (-l)'=(2fc-l)A! , (C30) 

Cl = l/k\ , (C31) 



(~<3 ^ ( _-\\ki+k2+l 2(fci + fc2) - 1 

'''''' ^ ' ki\{k2-2)\{ki+k2){ki+k2-l) 

+ (_l)fcl+fc2 + {-if A i + (-l)fci+fe2^(^l 



^1^2 A;i! (^2-2)! (A:i + A:2)(A:i + A:2-l) 



(C32) 
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Appendix D Check of the consistency of the mass— shell con- 



straint 



This Appendix proves the identity ( p.36 ). Let me exphcitly write out the operators 



>C, n, A and also 



Cfix,p) = p'^b^ f{x,p) + V— {F,,, f{x,p)} 

2 dpa 



(Dl) 



icfi 3 Ti^ ~ ~ 

Iif{x,p) = —p"^ [F^u, f{x,p)] - — D^D^ f{x,p) + 



4 dpo 



1. „ 1 



52 



1 



d 



dpadpp 4 ^ dp, 



(D2) 



Afix,p) 



ieti d 



8 dp. 



F,^,D^f{x,p) 



ieh ,. d"^ 
-p 



'■n 92 



8 " dp^dp^^^'-^'^'^^'^'P^^ 16 dp^dpp 



F^,Fl J{x,p) 



+ 



d^ 



Q3 



dpadpp 24 
d^ 



I,. 1. 



= P^Pa — fn ■ 



dpadppdp, 
d ' 



-+ 



f{x,p) 



(D3) 
(D4) 



To prove the identity ( p6D , one needs an expression for the commutator of two operators 
Da defined by Eq. ( 3.24| ). It can be shown that 



RaP^iuP'^ 



dPti 



where Vq, is the covariant and gauge invariant derivative operator. 
Therefore (compare with the formula (|C1| ) ), 



(D5) 



~ ~ 1 Z6 d 

F>a,Dp\^ f{x,p) = — [Faf3,f{x,p)] - Ra/SfiuP" f (x , p) 



(D6) 
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and 



P 



IP 7 P r ~ 



+ 



d 



+ f 2i?„^^,pV^^" - Rc 



Then, by noting the following property: 



{i, 







(D8) 



one can prove after elaborate calculations that the identity ( ^.361 ) is right up to the terms of 
the next adiabatic order. 
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